Ab-initio study of the vibrational properties of crystalline TeC^: the a, (3, and 7 phases 
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Based on density functional perturbation theory, we have studied the vibrational properties of 
three crystalline phases of tellurium dioxide: paratellurite «-Te02, tellurite /3-Te02 and the new 
phase 7-Te02, recently identified experimentally. Calculated Raman and IR spectra are in good 
agreement with available experimental data. The vibrational spectra of a- and /3-TeC>2 can be 
interpreted in terms of vibrations of Te02 molecular units. 



I. INTRODUCTION 

Tellurium oxide (TeC^) and TeCVbased glasses are 
promising active materials for optical switching devices 
due to their large non-linear polarizability 1 . and for opti- 
cal amplifiers due to their large cross section for Stim- 
ulated Raman Scattering (SRS). 2,3 Three crystalline 
phases of Te02 are well documented. The paratellu- 
rite a-TeO^ and the tellurite /J-TeO^ phases have been 
known for a long time. Recently, a third crystalline poly- 
morph metastable at normal conditions, 7-Te02r^ has 
been identified by x-ray powder diffraction of recrystal- 
lized tellurite glasses doped with metal oxides. The three 
phases can be described as different arrangements of cor- 
ner sharing TeC>4 unit (the distorted trigonal bipyramids 
in Fig. [T]) . These units have two longer Te — bonds 
(axial bonds) and two shorter ones (equatorial bonds) 
whose length varies in the different phases. Paratellu- 
rite, a-Te02, has a 3-dimensional fully connected struc- 
ture, while /3-Te02 has a layered structure with weakly 
bonded sheets. In 7~Te02, one Te — O bond is substan- 
tially larger than the other three and by breaking this 
latter bond in a 3-dimensional visualization of the net- 
work, a chain-like structure appears. The 7-Te02 phase 
has then been described in terms of a polymeric form 
made of TeC>3 unitsj&i By comparing the Raman spec- 
tra of TeC»2 glasses with those of the crystalline phases, 
it has also been proposed that the chain-like structure 
of 7~Te02 would represent the main structural feature 
of the glass4 Theoretical lattice dynamics calculations 
of the phonon spectra of crystalline TeC>2 have already 
appeared in lit er at ur e jL&S i 1 1 1 1 based on empirical in- 
teratomic force constants from two- and three body in- 
teractions or shell models^ The lattice dynamics cal- 
culations revealed that the phononic spectra of a-Te02 
and /3-Te02 can be well interpreted in terms of vibra- 
tions of weakly coupled TeC>2 molecules 1 ' 1 1 Accord- 
ingly, tellurite and paratellurite would be seen as molec- 
ular crystals made of TeC>2 molecules. The polymeric 
7-TCO2 would represent a chain-like polymerization of 
TeC>2 molecules giving rise to a structure more connected 
than paratellurite and tellurite phases. 

In this paper, we investigate further the three crys- 
talline phases of TeC>2 (a, f3, 7) by computing their struc- 
tural and vibrational properties from first principles. We 
aim at providing a compelling assignment of the exper- 
imental Raman and IR peaks to specific phonons which 



would allow to identify the vibrational signature of the 
different structural units in the crystals and in the glass. 

After a brief description of our theoretical framework 
in section II, we present our results on the structural and 
electronic properties in section III and on the vibrational 
properties in section IV. Section V reports our conclu- 
sions. 



II. COMPUTATIONAL DETAILS 



Calculations are performed within the framework of 
Density Functional Theory (DFT) within the simple 
local density approximation (LDA) and with gradient 
corrected exchange and correlation energy functionals 
(PBE 12 and BLYP 13 ), as implemented in the codes 
PWSCF and PHONONS H Calculations are performed 
with ultrasoft pseudopotentials^ Kohn-Sham (KS) or- 
bitals are expanded in a plane waves basis up to a kinetic 
cutoff of 30 Ry. Brillouin Zone (BZ) integration has been 
performed over Monkhorst-Pack (MP) meshes of 8x8x8, 
4x4x4 and 6x6x6 for a-, (3-, and 7-Te02, respectively*^ 
Equilibrium geometries have been obtained by optimiz- 
ing the internal and lattice structural parameters at sev- 
eral volumes and fitting the energy versus volume data 
with a Murnaghan function^ Residual anisotropy in the 
stress tensor at the optimized lattice parameters at each 
volume is below 0.6 kbar. Infrared and Raman spectra 
are obtained from effective charges, dielectric susceptibil- 
ities and phonons at the T point within density functional 
perturbation theory^ Effective charges have been com- 
puted with norm-conserving pseudopotentials, 1 ^ and a 
kinetic cutoff of 70 Ry. Relevant formula for the calcula- 
tion of the IR and Raman spectra are given in section IV. 
Preliminary tests on the transferability of the pseudopo- 
tcntials and on the accuracy of the exchange and corre- 
lation functionals have been performed on the structure 
and vibrational frequencies of the Te02 molecule. Re- 
sults are presented in Table Q] for PBE and BLYP func- 
tionals and compared with experimental data and previ- 
ous B3LYP all-electrons calculations. 
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FIG. 1: The Te04 bipyramidal unit, building block of the 
three TeC>2 polymorphs. The two small circles indicate the 
Te lone-pair. Dark (light) spheres indicate Te (O) atoms. 




TABLE II: Structural parameters of a— Te02, calculated 
with LDA functional, with PBE functional at the theoretical 
(PBE) and experimental (PBEexp) lattice parameters, and 
with BLYP functional at the experimental lattice parameters 
(BLYPexp). Experimental X-ray data (Exp) are from Ref.— . 





LDA 


PBE 


PBEexp 


BLYPexp 


Exp 


Cell paramet 


ers (A) 










a 


4.819 


4.990 






4.8082 


c 


7.338 


7.546 






7.612 


Atomic posit 


ions 










Te x 


0.0025 


0.0272 


0.0156 


0.0201 


0.0268 


O a 


0.1402 


0.1467 


0.1371 


0.1388 


0.1368 


o p 


0.2379 


0.2482 


0.2492 


0.2519 


0.2576 


O 7 


0.2065 


0.1968 


0.1954 


0.1926 


0.1862 



TABLE I: Bond length, bond angle and vibrational frequen- 
cies of Te02 molecule, compared with previous B3LYP>2i and 
experimental dataS. Bi and Ai(2) are symmetric and an- 
tisymmetric stretching modes, respectively, while Ai(l) is a 
bending mode of the TeC>2 molecule. 









0J 


(cm" 1 ) 






Te-0 (A) 


Te-O-Te (deg.) 






Si 


PBE 


1.832 


111.7 


266 


804 


806 


BLYP 


1.820 


111.8 


265 


813 


824 


B3LYP 


1.844 


112.9 


250 


883 


921 


Exp 


1.840 


112 


294 


810 


849 



III. STRUCTURAL PROPERTIES 
A. Q-Te0 2 

The paratellurite crystal (space group P4i2i2, Z)|, n. 
92) has tetragonal symmetry with four formula units per 
unit cell. Two atoms are independent by symmetry: one 
Te atom at position (x, x, 0) and one oxygen atom at po- 
sition (a, (3, 7). The experimental and theoretical (LDA 
and PBE) lattice parameters and position of the symme- 
try independent atoms are given in table [Til As usual the 
LDA (GGA) underestimates (overestimates) the equilib- 
rium volume with respect to experimental data. The 
bulk modulus and its derivative with respect to pres- 
sure obtained from the fitted Murnaghan equation of 
state (PBE functional) are B=33.0 GPa and B'=8.1 
to be compared with the experimental value of i?=45.3 
GPa. 20 These results are in good agreement with previ- 
ous all-electron calculations which confirms the transfer- 
ability of our pseudopotentials. 2 i Concerning the inter- 
nal structure of a-Te02, a Te atom is coordinated with 
four oxygen atoms, two equatorial atoms at closer dis- 
tance (exp. 1.88 A) and two axial atoms at larger dis- 
tance (exp. 2.12 A) 4. The structure can be seen as a 
network of corner-sharing Te04 units (Fig. [SJ) , the oxy- 
gen atoms bridging two Te atoms. The network is formed 



by six-membered rings of Te atoms with bridging oxygen 
atoms. As already mentioned, lattice dynamics calcu- 
lations with empirical interatomic potentials suggest an 
alternative picture of a-Te02 as formed by weakly inter- 
acting Te02 molecules i 8 ' 10 ! 11 The calculated independent 
bond lengths and angles are compared with experimen- 
tal data in Table IIIII for LDA and PBE functionals at 
the theoretical equilibrium lattice parameters, and for 
PBE and BLYP functionals at the experimental lattice 
parameters. Although the error in the lattice parameters 
is within the expected accuracy for DFT calculations, the 
error in the bond lengths at the theoretical lattice param- 
eters is somehow larger than usual (cfr. Table HIT)) . 

In particular the calculated equatorial and axial 
Te — O bond lengths are much closer to each other than 
in the experimental structure. This misfit is larger for 
LDA than for PBE functional, but it is sizably reduced 
by optimizing the internal geometry at the experimental 
lattice parameters with a marginal improvement for the 
BLYP functional (cfr. Table [TTI]). The change in the lat- 
tice parameters, bond lengths and bond angles with vol- 
ume is reported in Fig. [3] The increase of the c/a ratio 
with pressure is in agreement with the experimental data 
of Refill. By increasing the volume, the length of the 
shorter equatorial bonds decreases while the longer axial 
bonds increases in length. This behavior would support 
the picture of a-Te02 as made of Te02 molecules, in- 
ferred from previous lattice dynamics calculation a 8 ' 10 ! 11 
which turns out to be in good agreement with our ab- 
initio phonon calculations reported in Sec. IV. Although 
useful for the analysis of the vibrational spectra, the pic- 
ture of paratellurite as a molecular crystal does not ac- 
count for all the properties of a-Te02- For instance, it 
can be hardly conciliated with the mostly ionic character 
and electronic bandwidths of several eV emerged from 
previous ab-initio calculations^^. 
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FIG. 2: Structure of a-TeC>2. (a) Perspective view, two unit 
cells are shown; 02 and 03 (Ol and 04) label the equatorial 
(axial) oxygen atoms, (b) Four unit cells are drawn along c 




FIG. 3: (a) Lattice parameters (a, c), (b) bond lengths and 
(c) bond angles of Q-Te02 as a function of unit cell volume. 
The arrow indicates the theoretical equilibrium value. The 
series in panel (a) correspond (from bottom to top) to a and 
c axes. Te — 03 and Te — 01 are the bond lengths of the axial 
and equatorial bonds, respectively. 




Volume [A ValJ 




B. /?-Te0 2 

The tellurite crystal (space group Pbca, -D|5, n. 61) 
has orthorhombic symmetry with eight formula units per 
unit cell.— Two oxygen and one tellurium atoms are inde- 
pendent by symmetry. The experimental and theoretical 
(PBE) lattice parameters and positions of the symmetry 
independent atoms are given in table [IV] The struc- 
tural unit of /3-Te02 is still a Te04 trigonal bipyramid, 
as in ev-Te02, which forms a network by sharing ver- 
tices (bridging oxygen atoms). However, as opposed to 
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TABLE III: Bond lengths and angles for a— TeC>2, calculated 
with LDA functional, with PBE functional at the theoretical 
(PBE) and experimental (PBEexp) lattice parameters, and 
with BLYP functional at the experimental lattice parameters 
(BLYPexp). Experimental X-ray data (Exp) are from Ref.—. 
Atoms are labelled according to Fig. [2] 





LDA 


PBE 


PBEexp 


BLYPexp 


Exp 


Bond lengths (A) 


Te-Ol 


2.028 


1.944 


1.954 


1.924 


1.878 


Te-03 


2.170 


2.118 


2.145 


2.134 


2.122 


Te-Te' 


3.855 


3.838 


3.806 


3.779 


3.742 


Angles (dej 


;rees) 










Ol-Te-02 


99.8 


103.6 


103.6 


103.7 


103.4 


03-Te-04 


161.5 


171.2 


165.3 


166.7 


168.0 


Te-02-Te' 


135.0 


137.1 


136.4 


137.2 


138.6 



TABLE IV: Structural parameters of f3— Te02, computed us- 
ing PBE functional at the theoretical (PBE) and experimental 
(PBEexp) lattice parameters, and with BLYP functional at 
the experimental lattice parameters (BLYPexp). Experimen- 
tal X-ray data (Exp) are from Ref.—. 





PBE 


PBEexp 


BLYPexp 


Exp 


Cell parameters 


(A) 








a 


12.232 






12.035 


b 


5.523 






5.464 


c 


5.803 






5.607 


Atomic positions 


Te x 


0.1173 


0.1193 


0.1189 


0.1181 


Tey 


0.0263 


0.0184 


0.0211 


0.0252 


Te z 


0.3748 


0.3677 


0.3682 


0.3378 


0(1) X 


0.0314 


0.3235 


0.0318 


0.028 


O(l) y 


0.6491 


0.6422 


0.6412 


0.634 


0(1) z 


0.1613 


0.1721 


0.1703 


0.171 


0(2) x 


0.1714 


0.1743 


0.1723 


0.168 


0(2) y 


0.2174 


0.2108 


0.2124 


0.221 


0(2) z 


0.0812 


0.0694 


0.0714 


0.081 



paratellurite, the /9-Te02 crystal has a layered structure 
with layers oriented perpendicular to the a axis (Fig. [4|) 
and an interlayer distance of ~ 3 AJi We can still rec- 
ognize two longer axial bonds and two shorter equatorial 
bonds, but here all the four Te — O bonds of the TeC>4 
unit differ in length. The 2D network is formed by two- 
and six-membered Te rings. Still, /3-Te02 as well can 
be seen as a molecular crystal made of TeC>2 molecules 
according to lattice dynamics calculations 10 ' 1 1 

The calculated independent bond lengths and angles 
are compared with experimental data in Table fVl for PBE 
functional at the theoretical equilibrium lattice parame- 
ters, and for PBE and BLYP functionals at the experi- 
mental lattice parameters. The same comments on the 
internal structure presented for a-Te02 holds here for /3- 



FIG. 4: Layered structure of /3-Te02. (a) A projection along 
ab plane, (b) Sketch of a single layer. Ol and 02 (Ol' and 
02') label equatorial (axial) atoms. 




TeC>2. The spread in the bond lengths of the theoretical 
geometry at the theoretical lattice parameters (PBE) is 
reduced with respect to experiments. The misfit is re- 
duced by optimizing the internal geometry at the experi- 
mental lattice parameters, with a marginal improvement 
for the BLYP over PBE functional. The error in the 
internal structure thus comes primarily from the overes- 
timation of the equilibrium volume (6 %) by PBE func- 
tional. A possible source of error might be the neglect of 
van der Waals interactions among highly polarizable Te 
atoms, absent in LDA and GGA functionals, whose in- 
clusion might bring theoretical and experimental lattice 
parameters to a closer agreement. The bulk modulus and 
its derivative with respect to pressure obtained from the 
fitted Murnaghan equation of state (PBE functional) are 
5=18.4 GPa (B=33 GPa for a-Te0 2 ) and B'=18.1. The 
larger compressibility of /3-Te02 with respect to a-Te02 
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TABLE V: Bond lengths and angles for /3-Te0 2 , using PBE 
functional at the theoretical (PBE) and experimental (PBE- 
exp) lattice parameters, and with BLYP functional at the 
experimental lattice parameters (BLYPexp). Experimental 
X-ray data (Exp) are from Ref.— . Atoms are labelled accord- 
ing to Fig. |1 





PBE 


PBEexp 


BLYPexp 


Exp 


Bond lengths (A) 


Te-Ol 


1.953 


1.960 


1.941 


1.876 


Te-02 


1.968 


1.975 


1.967 


1.893 


Te-Ol' 


2.192 


2.186 


2.181 


2.153 


Te-02' 


2.110 


2.083 


2.067 


2.068 


Te-Te' 


3.232 


3.240 


3.230 


3.168 


Angles (degrees) 


01- Te-02 


97.4 


96.8 


96.9 


98.8 


01'-Te-02' 


169.5 


168.9 


168.6 


165.9 


Te-Ol-Te' 


102.3 


102.6 


103.0 


103.5 



is due to the weak interlayer interaction, as demonstrated 
by the evolution of the lattice parameters with volume 
reported in Fig. [5^,. The a axes, perpendicular to the 
layers, undergoes larger changes with pressure than the 
other two axes do. However the overestimation of the lat- 
tice parameters by the PBE functional is almost the same 
for all the three axes which implies that the interlayer in- 
teraction should have a strong electrostatic contribution. 
Would the interlayer cohesion be simply given by van der 
Waals interactions, absent in GGA functional, the over- 
estimation of the a axis by the PBE functional should 
have been much larger. The change of bond lengths and 
bond angles with volume is reported in Fig. [5p-c. By 
increasing the volume, the length of the shorter equato- 
rial bonds decreases while the longer axial bonds increase 
in length. This behavior would support the picture of /?- 
Te02 as made of TeC>2 molecules, analogously to a-Te02- 



C. 7-Te0 2 

This phase has been recently identified by recrystalliz- 
ing amorphous TeC>2 doped with oxide (Nb, W oxides, 10 
% molar) The structure is orthorhombic (space group 
P2i2i2i, D2, n. 19) with four formula units and three 
independent atoms per cell. The theoretical (PBE) and 
experimental lattice parameters and the positions of the 
independent atoms are reported in table IVIl The error 
in the theoretical equilibrium volume (11 %) and lattice 
parameters (3 %) is slightly larger for 7-Te02 than for 
the other phases. The structure of 7-Te02 can still be 
seen as formed by corner sharing TeC>4 units with two- 
fold coordinated oxygen atoms and four-fold coordinated 
Te atoms. However, the difference in length of the four 
Te — O bonds is 0.34 A, a value sizably larger than in 
the other phases (0.27 A for a-Te0 2 and 0.24 A for 
/3-Te0 2 ). By considering the longest Te - O bond (2.198 



FIG. 5: (a) Lattice parameters (a, b, c), (b) bond lengths and 
(c) bond angles of /3-Te02 as a function of unit cell volume. 
The arrow indicates the theoretical equilibrium value. The 
three series in panel (a) correspond (from bottom to top) to 
b, c and a axes. The atom labels refer to Fig. U Te — Ol and 
Te - 02 are equatorial bonds, while Te - Ol' and Te - 02' 
are axial bonds. 
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TABLE VI: Structural parameters of 7-Te0 2 , using PBE 
functional at the theoretical (PBE) and experimental (PBE- 
exp) lattice parameters, and with BLYP functional at the 
experimental lattice parameters (BLYPexp). Experimental 
X-ray data (Exp) are from Ref.— . 





PBE 


PBEexp 


BLYPexp 


Exp 


Cell parameters 


(A) 








a 


5.176 






4.898 


b 


8.797 






8.576 


c 


4.467 






4.351 


Atomic positions 


Te x 


0.9581 


0.9797 


0.9771 


0.9696 


Tey 


0.1032 


0.1003 


0.1007 


0.1016 


Te z 


0.1184 


0.1407 


0.1380 


0.1358 


0(1) X 


0.7641 


0.7827 


0.7801 


0.759 


0(1) y 


0.2851 


0.2905 


0.2888 


0.281 


0(1) z 


0.1645 


0.1853 


0.1805 


0.173 


0(2) x 


0.8599 


0.8656 


0.8646 


0.855 


0(2) y 


0.0406 


0.0371 


0.0361 


0.036 


0(2) z 


0.7131 


0.7237 


0.7278 


0.727 



A) as a non-bonding interaction, 7-Te02 would appear 
as a network of Te03 units with Te atoms three-fold co- 
ordinated, two oxygen atoms two-fold coordinated and 
one non-bridging one-fold coordinated oxygen atom. 

The structure of 7-TeC>2 depicted as a fully connected 
3D network or as a one-dimensional polymeric structure 
of TeC>3 units is shown in Fig. [6^l and Fig. [BJd, respec- 
tively. The calculated independent bond lengths and an- 
gles are compared with experimental data in Table IVIII 
for PBE functional at the theoretical equilibrium lattice 
parameters and for PBE and BLYP functionals at the 
experimental lattice parameters. Still, the spread in the 
bond lengths is reduced in the theoretical geometry at 
the theoretical lattice parameters (PBE) with respect to 
experiments. The larger error in the theoretical lattice 
parameters shows up in larger errors in the Te — O bond 
lengths in 7-TeC>2 with respect to the other phases (cfr. 
tables |nT] and [Vj) . The misfit is largely reduced by opti- 
mizing the geometry at the experimental lattice param- 
eters, with a marginal improvement for the BLYP over 
PBE functional. In particular, the optimized structure 
at the experimental lattice parameters reproduces the 
spread of 0.34 A in the Te — O bond lengths. The change 
in the lattice parameters, bond lengths and bond angles 
with volume is reported in Fig. [7J The shorter Te — 01 
bond does not vary with pressure, while both the equa- 
torial Te — 02 bond and the longer Te — 01' axial bond 
decrease in length by increasing the volume. Regard- 
ing the picture of 7-Te02 as a polymeric phase, we note 
that the longest Te - O bond of the Te0 4 unit (2.198 
A) is only 5 % longer than the other axial bond. On 
the other hand, one notes that in 7~Te02 the difference 
in length between the longer equatorial bond and the 



FIG. 6: Structure of 7-Te02- Two unit cells are drawn, 
(a) Three-dimensional structure, (b) Chain-like structure ob- 
tained by breaking the longest Te — O bond of the Te04 units 
of panel (a). 01 and 02 (01' and 02') label equatorial (axial) 
atoms. 




short axial bond is the smallest among the three crys- 
talline phases. The polymeric structure would then be 
seen as generated from the formation of strong inter- 
molecular bonds between the Te02 structural units of 
tellurite and paratellurite rather than from the breaking 
of one bond of the Te04 units^i^ The one-dimensional 
character of 7-Te02 has to be supported by the observa- 
tion of strong anisotropics in the physical properties of 
this phase. To address this issue, we have calculated the 
elastic constants of 7-Te02 from the difference in the to- 
tal energy of strained crystals. A maximum strain of 3 % 
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TABLE VII: Bond lengths and angles for 7-TeCb, using PBE 
functional at the theoretical (PBE) and experimental (PBE- 
exp) lattice parameters, and with BLYP functional at the 
experimental lattice parameters (BLYPexp). Experimental 
X-ray data (Exp) are from Ref.— . Atoms are labelled accord- 
ing to Fig. [fj] 





PBE 


PBEexp 


BLYPexp 


Exp 


Bond lengths (A) 


Te-Ol 


1.900 


1.905 


1.889 


1.859 


Te-02 


1.960 


1.974 


1.949 


1.949 


Te-Ol' 


2.252 


2.257 


2.240 


2.198 


Te-02' 


2.119 


2.092 


2.081 


2.019 


Te-Te' 


3.596 


3.571 


3.559 


3.521 


Angles (degrees) 


01- Te-02 


99.1 


100.7 


100.9 


101.5 


01'-Te-02' 


156.8 


150.9 


152.5 


153.6 


Te-Ol-Te' 


123.6 


122.8 


124.0 


125.0 



has been considered. For stretching along the principal 
axes we obtain Cn—38 GPa, C*22=41 GPa, and 6*33=43 
GPa, where the Voight notation has been used. The bulk 
modulus and its derivative with respect to pressure ob- 
tained from the fitted Murnaghan equation of state (PBE 
functional) are £=16 GPa (B=33 GPa for a-Te0 2 ) and 
B'=6.2. The elastic properties are thus nearly isotropic 
and do not support the picture of 7~Te02 as made of one- 
dimensional polymeric structure. The clastic constants 
reported above differ from the results of calculations with 
empirical interatomic potentials (C*n=30 GPa, 6*22=40 
GPa, and 6*33=67 GPa, Refi), mainly in the C33 com- 
ponent which shows an appreciable lower compressibility 
along the chain (z) axis, not confirmed by our ab-initio 
results. 

We have also contrasted the electronic band structure 
of 7-Te02 and /3-Te02 phases in Fig. [S] The electronic 
band structure of a-Te02 has been computed within 
DFT in a previous work.— In /3-Te02, the band disper- 
sion along the direction T-X, perpendicular to the 2D 
layers, is sizably smaller than along the T-Y and T-Z di- 
rections within the layer plane. Conversely, in 7-Te02, 
the dispersion is similar along the three orthogonal di- 
rections T-X, T-Y, T-Z. Nor the dispersion along T-Z is 
much larger in 7~Te02 than in /3-Te02, as one would ex- 
pect by comparing a polymeric phase with chains along 
z (7 phase) with a mainly molecular crystal made of 
Te02 molecules ((3 phase). Therefore, the electronic band 
structure does not provide evidence of a polymeric nature 
of7-Te0 2 . 

We have also computed formation enthalpy and forma- 
tion free energy at normal conditions (298.15 K, p=l bar) 
for the three cristalline phases, by including the phononic 
contribution to free energy and enthalpy. Only T-point 
phonons have been included in the calculation (see sec- 
tion UVj. For a-Te02 we have verified that the phononic 
contribution to the free energy of -7.8 kJ/mol, which in- 



FIG. 7: (a) Lattice parameters (a, b, c), (b) bond lengths and 
(c) bond angles of 7-Te02 as a function of unit cell volume. 
The arrow indicates the theoretical equilibrium value. The 
three series in panel (a) correspond (from bottom to top) to 
c, a and b axes. The atom labels refer to Fig. HJ] 
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FIG. 8: Electronic band structure of (a) /3-Te02 and (b) 
7-Te02. The orthorhombic irreducible Brillouin zone is re- 
ported in the inset, equal in shape (but not in size) for both 
crystals. 



FIG. 9: Angular dispersion of E and A2 IR-active modes in 
a-Te02- Dashed lines represent unshifted (i.e. TO) E modes. 
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eludes only T-point phonons, changes to -6.7 kJ/mol by 
including in the calculation the full phonon density of 
states. Formation free energies and enthalpies have been 
computed with PBE functionals at the theoretical equi- 
librium structure for all the three phases. 

The calculated formation enthalpy (AHf) and free 
energy (AG/) of a-Te02, with respect to crystalline 
Te and gaseous O2 (the standard states), are AHf =- 
287.5 kJ/mol (exp. -322.6 kJ/mol) 22 , and AG/ =-232.8 
kJ/mol (exp. -270.3 kJ/mol)22. The error is within 
the usual accuracy of DFT methods. The free energies 
of gaseous Te2 and O2 have been calculated from the 
molecule optimized in the triplet stated by including 
translational, vibrational and rotational contributions. 
The free energy of crystalline Te has been computed 
from theoretical and experimental 2 ^ free energies of the 
gaseous Te2 molecule. 

For /3-Te02 the calculated formation enthalpy and 
free energy are AHf =-289.3 kJ/mol and AG/ =-236.3 
kJ/mol, while for 7-Te0 2 AH f =-284.6 kJ/mol and 
AG/ =-229.8 kJ/mol. To our knowledge, no experimen- 
tal data on formation energies are available for /3-Te02 
and 7-Te02- 
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IV. VIBRATIONAL PROPERTIES 

Phonons at the T-point have been calculated within 
density functional perturbation theory^ The modes 
which display a dipole moment couple to the inner macro- 
scopic longitudinal electric field which shifts the LO 
phonon frequencies via the non-analytic contribution to 
the dynamical matrix 18 



4tt Z aa , (n)q al Zpp, {n^qp, 

v q ■ e°° • q 



(i) 



where Z_ and e°° are the effective charges and electronic 
dielectric tensors, V is the unit cell volume and q is the 
phononic wavevector. The macroscopic field contribution 
to the dynamical matrix (Eq. [T]) introduces an angular 
dispersion of the phonons at the T point, i.e. the limit 
of the phononic bands uj(q) for q — > depends on the 
angles formed by q with the principal axis. For a uniax- 
ial crystal as a-Te02, the phonon frequencies depend on 
the angle 8 formed by q with the optical axis. For the 
orthorhombic crystals /3-Te02 and 7-Te02, the phonon 
frequencies depend on the polar angles 9 and <p which 
define the orientation of q with respect to the crystallo- 
graphic axes, coinciding with the optical principal axes. 

The dielectric tensor is given in terms of phonons and 
effective charges by 



4-7T 3N 1 

e af3 {u)) = + TT S~~] Pa (j)Pfl (j) 2 

v — ' ~~ 



where p(j) is given by 



A' 



(2) 



(3) 



£=1 



The absorption coefficient a 7 for light polarized along 
the crystallographic axis 7 is 



solid angle of the possible wave vectors q of the trans- 
mitted wave. The absorption coefficient for a generic q 
is obtained by solving the Fresnel's equations for the di- 
electric tensor defined by Eq. 

The differential cross section for Raman scattering 
(Stokes) in non-resonant conditions is given by the fol- 
lowing expression (for a unit volume of scattering sample) 
4 

= ^ \ e s ■ B? ■ e L \ 2 (n fl (w) + l)*(w - u)j), 



dfldui 



j 



(5) 

where ub(uj) is the Bose factor, uis is the frequency of 
the scattered light, e$ and are the polarization vec- 
tors of the scattered and incident light, respectively 27 ^ 8 
The Raman tensor R? associated with the j-th phonon 
is given by 



R 3 




(6) 



where V a is the unit cell volume, r(/c) is the position of 
the K-th atom and x°° — (i 00 ~ t j/Air is the electronic 

susceptibility. The tensor R J is computed from x°° by 
finite differences, by moving the atoms independent by 
symmetry with maximum displacement of 0.0 f A. 

Whenever the experimental Raman spectrum is avail- 
able only for a polycrystalline sample, Eq. [5] must be 
integrated over the solid angle in order to compare the 
theoretical spectra with experiments. In particular, the 
Raman spectrum for non-polarized light is obtained by 
summing over all possible polarization vectors e$ and 
consistent with the scattering geometry. In the case 
of /3-Te02, the inner longitudinal macroscopic field has 
no effect on the Raman active modes due to the pres- 
ence of an inversion symmetry. In fact, the Raman ac- 
tive g-modes have no dipole moment to couple with the 
macroscopic longitudinal field. As a consequence for /?- 
TeC>2 there is no angular dispersion for the Raman-active 
modes and the total cross section for unpolarizcd light in 
backscattcring geometry is obtained from Eq. [5] with the 
substitution 



a~,(uS) = — Im Cy-yiu) + in, rj — > 0) 
nc 



(4) 



where c is the velocity of light in vacuum and n(u>) 
is the (frequency dependent) real part of the refractive 
index. 

In order to compare the theoretical results with avail- 
able experimental data we have also computed the IR 
absorption for a polycrystalline sample. This amounts to 
average the absorption coefficient for the two polariza- 
tions of the electromagnetic wave (ordinary and extraor- 
dinary waves for uniaxial crystals, for instance) over the 



4(-^xa; 4~ -^"yy 4~ ^zz) 4~ ^{^xy ^xz 4~ ^yz 



{BxxRyy ~\- RxxRz 



RzzRyy) 



30\e s -R j -e L | 2 (7) 



Conversely, for a-Te02 and 7~Te02 the presence of an- 
gular dispersion requires the integral over the solid angle 
to be performed by summing over discrete points in the 
polar angles. Phonon frequencies and eigenvectors enter- 
ing in Eq. [5] and [H] are modified by the macroscopic lon- 
gitudinal field along the direction assigned by the crystal 
momentum transferred in the scattering process. 

The (^-functions in Eq. 2] and [5] are approximated by 
Lorentzian functions as 



5(w - Luj) = - 



7T (W 2 



,2\2 



(8) 
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with rjj fitted on the experimental Raman peaks when 
available and otherwise assigned to a constant value of 4 
cm" 1 (cfr. Table 



a-Te0 2 



Phonons at the T-point can be classified according to 
the irreducible representations of the -D4 point group of 
a-Te0 2 as T = 4Ai + AA 2 + 5Bi + AB 2 + 8E, where 
the acoustic modes (A 2 and E) have been omitted. The 
Ai, Bx, B 2 and E modes and Raman active, while the 
A 2 and E modes are IR active. The calculated phonon 
frequencies at the L point neglecting the contribution of 
the macroscopic longitudinal field are given in Table lyTTTl 
for the BLYP calculation at the experimental lattice pa- 
rameters. The angular dispersion of A 2 and E modes 
due to the macroscopic electric field is reported in Fig. 
IH1 it is in good agreement with the angular dispersion 
obtained experimentally from the IR absorption peaks of 
the extraordinary waves (cfr. Fig. 3 of Ref.— ). Experi- 
mental frequencies and IR activities of IR-active modes 
(A 2 and E(TO) modes) are also reported in Table IVIlTl 
Experimental frequencies from Raman spectra are also 
reported in Table IVlTll for A\, B\ and B 2 modes. For 
Raman active E(LO) modes the comparison is made di- 
rectly between the theoretical and experimental Raman 
spectra, since their frequencies depend on the scatter- 
ing geometries. Phonon frequencies calculated with PBE 
or LDA functional at the theoretical lattice parameters 
and with PBE and BLYP functionals at the experimental 
lattice parameters are compared in Table rvTnl The best 
agreement with experiments is obtained for BLYP calcu- 
lations at the experimental lattice parameters which give 
the equilibrium internal geometry closest to the exper- 
imental structure (cfr. section IIIA). The IR oscillator 
strengths are in good agreement with experiments, apart 
from the overestimation of the intensity of mode A 2 (2). 
Experimental data from a more recent work 30 are very 
similar (within 5 cm -1 ) to the older experimental data 
reported in Table IVlTTI but for the oscillator strength (/) 
of mode A 2 {2) which is 3.33 in Ref^, in closer agreement 
with our theoretical result. The IR absorption spectrum 
for a polycrystalline sample, computed as described in 
the previous section is reported in Fig. [121 it nas to be 
compared with the experimental spectrum reported in 
Fig. 2 of Ref. 8 . The linewidth of the phononic modes 
is chosen according to the IR experimental data on sin- 
gle crystal^ The two broad bands and the shoulder at 
lower frequency in the experimental spectrum^ are well 
reproduced, but for the already mentioned redshift in the 
theoretical frequencies. 

The Raman tensor (Eq. [5]) for the Raman-active irre- 
ducible representations has the following form^i 



TABLE VIII: Theoretical phonon frequencies of «-TeC>2 at 
the r point, oscillator strengths (fj in Eq. Q of IR active 
modes and squared coefficients of the Raman tensor of the 
Raman active modes, a 2 and b 2 for Ai, c 2 for Bx, d 2 for B 2 , 
and e 2 for E x , E y modes (in units of 10 _3 A 3 , see section IV). 
The contribution of the inner longitudinal macroscopic field 
is not included (LO-TO splitting). All the theoretical val- 
ues correspond to calculations with the BLYP functional at 
the experimental lattice parameters. Experimental data on 
phonon frequencies and oscillator strengths from IR absorp- 
tion spectra at 85 K— are given in parenthesis for (TO) A2 
and E modes. Experimental phonon frequencies of Ax, Bx 
and B2 modes (in parenthesis) are taken from Raman spectra 
at 85 K 33 . Alo indicate the phonon linewidths obtained by 
fitting the main peaks of the polycrystalline Raman spectra 
at 295 K (cfr. Fig. 



Mode 


ui (cm 1 ) 


Si 


a (c , a , e ) 


b 




Bi (1) 


42 


(62) 




2.044 






A ( 1 \ 

Ai (1) 


90 


(76) 7.758 (12.95) 








E (1) 


124 


(124) 0.919 


(1.28) 


5.158 




2.9 


Bi (2) 


128 






0.042 






A X (1) 


143 


(152) 




4.939 


7.561 


5.8 


B 2 (1) 


148 


(157) 




2.767 




2.2 


Bx (3) 


175 


(179) 




0.324 






E(2) 


175 


(177) 5.667 


(8.05) 


0.980 






E(3) 


212 


(212) 3.127 


(2.02) 


0.043 






Ax (2) 


214 


(218) 




0.338 


0.058 




Bx (4) 


229 


(235) 




0.514 






A 2 (2) 


269 


(265) 3.156 


(0.77) 








B 2 (2) 


272 


(281) 




0.330 






E{A) 


291 


(299) 5.121 


(3.95) 


0.045 






A 2 (3) 


302 


(325) 4.096 


(4.70) 








E(5) 


325 


(336) 0.005 


(0.25) 


0.077 






E(6) 


336 


(379) 0.006 


(0.01) 


0.863 




6.4 


Ax (3) 


383 


(391) 




2.421 


7.923 


12.1 


B 2 (3) 


406 


(415) 




0.136 






A 2 (4) 


535 


(570) 3.802 


(3.86) 








Bx (5) 


537 


(589) 




1.551 




14.0 


E(7) 


594 


(644) 1.268 


(1.27) 


0.200 






Ax (4) 


598 


(649) 




10.242 


14.815 


4.7 


E(8) 


700 


(774) 0.225 


(0.21) 


0.671 






B 2 (4) 


726 


(786) 




0.574 







Ax 



B> 



a 

. b 



Bx 



(-e) 
e 



(-e) e 
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The coefficients a,6,c,d, and e calculated from first 
principles as outlined above are given for each mode 
in Table I Villi The calculated derivatives of the di- 
electric tensor with respect to the displacement of the 
atoms independent by symmetry are given as additional 
materials.— 

The theoretical Raman spectra for different scatter- 
ing conditions, each selecting particular modes, are com- 
pared with experimental data^^ in Fig. [101 The the- 
oretical spectra include the shift in phonon frequencies 
due to the contribution of the macroscopic longitudinal 
field along the direction assigned by the crystal momen- 
tum transferred in the scattering geometry. The Raman 
spectrum for a polycrystalline sample is also computed 
and compared with experimental data in Fig. 1131 For 
each wavevector q transfered in the scattering process 
the phonon eigenvector and frequency is computed by 
including the non-analytic part of the dynamical matrix. 
The solid angle average is performed over 256 angles in- 
dependent by symmetry. The linewidth of the most in- 
tense Raman peaks in the theoretical spectra has been 
obtained by fitting the experimental spectra with a sum 
of Lorentzian functions, either at 85 K for the single crys- 
tal spectra (Fig. [101) or a * 295 K for the polycrystalline 
spectra. The high temperature linewidth are given in Ta- 
ble I Villi The linewidth of the weaker or less identifiable 
peaks is set to 4 cm -1 . 

The calculated Raman spectra are in good agreement 
with experiments on peak intensities and also on peak 
frequencies in the low-frequency region, while the high- 
frequency modes appear too soft by nearly 8% (see Table 
IVIII[) . The discrepancy can possibly be ascribed to the 
partial homogenization of axial and equatorial bonds re- 
sulting from the use of approximate exchange-correlation 
functionals. This effect is smaller for BLYP than for PBE 
functionals (cfr. section IIIA). The experimental Raman 
peak at 218 cm -1 assigned to a B\ mode in Ref,— can 
be better assigned to the ^4i(2) mode on the base of our 
calculations (cfr. Table |VTTT| which is also not inconsis- 
tent with experimental two-phonons IR spectra of Ref^. 
More on the displacement pattern of the different modes 
is given in section IVD where the Raman spectra of the 
three phases are compared. 

Phonon dispersion relations along the high symmetry 
directions of the Irreducible Brillouin Zone (IBZ) and 
phonon Density Of States (DOS) are reported in Fig. 
1111 The dynamical matrices have been computed for 30 
q points along the T — Z — R — X — T path shown in 
Fig. QTJ or on a 4x4x4 grid in q-space for the DOS 
calculation. A Fourier interpolation technique provided 
the dynamical matrix at the other points of the BZ. 18 
For the DOS calculation, we have used the tetrahedron 
method on a mesh of 20 x 20 x 20 points in the IBZ. 

An energy gap separates the low-frequency phonon 
branches from the eight highest frequency bands, which 
can be assigned to stretching modes of the eight inequiv- 
alent Te — O bonds per cell, i.e. symmetric and anti- 
symmetric stretching modes of the four Te02 in the unit 



FIG. 10: Raman spectra of «-TeC>2 single crystal for different 
scattering geometries. Theoretical spectra (th.) are com- 
puted with the BLYP functional at experimental lattice pa- 
rameters. The experimental spectra at 85 K are adapted from 
Ref.—. The character of the modes active in each scattering 
geometry is also reported. The linewidths of the theoretical 
peaks are fitted on the experimental spectra whenever the 
correspondence is unambiguous. The Raman intensities are 
reported on a logarithmic scale for sake of comparison with 
the experimental data of Ref.—. The experimental peak is 
normalized to the strongest peak of the corresponding theo- 
retical spectrum. 
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FIG. 11: Phonon dispersion relations and density of states of 
a-Te02, calculated with BLYP functional at the experimental 
lattice parameters. 

™ 1 i • 1 1 




FIG. 12: IR absorption (a(ui), cfr. Eq. [4} spectra for poly- 
crystalline samples of a-, /3-, and 7-Te02; linewidths in a- 
Te02 are taken from experimental data on single crystals 29 . 
In the lack of experimental data, the linewidths /3-, and 7- 
Te02 are chosen analogously to those of a phase, i.e. 5 cm -1 
for modes with frequency lower than 150 cm' 1 , and 15 cm -1 
for higher frequency modes. 
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TABLE IX: Phonon frequencies (cm -1 ) of a-Te02 computed 
with BLYP functional at the experimental lattice parameters 
(BLYPexp), with PBE iand LDA functionals at the theoret- 
ical lattice parameters and with PBE at the experimental 
lattice parameters (PBEexp). 



Mode 


BLYPexp 


PBEexp 


PBE 


LDA 


Bi (1) 


42 


48 


62 


69 


A 2 (1) 


90 


94 


90 


106 


E (1) 


124 


120 


116 


142 


Bi (2) 


128 


127 


125 


129 


Ai (1) 


143 


140 


134 


135 


B 2 (1) 


148 


141 


136 


131 


Bi (3) 


175 


171 


150 


166 


E (2) 


175 


171 


162 


179 


E (3) 


212 


214 


193 


213 


Ai (2) 


214 


207 


195 


206 


Bi (4) 


229 


234 


218 


251 


A 2 (2) 


269 


267 


254 


256 


B 2 (2) 


272 


261 


235 


241 


E (4) 


291 


282 


262 


262 


A 2 (3) 


302 


292 


286 


284 


E (b) 




OlD 




OZ / 


E(6) 


336 


335 


317 


334 


Ai (3) 


383 


377 


364 


402 


B 2 (3) 


406 


399 


385 


404 


A 2 (4) 


535 


507 


515 


471 


Bi (5) 


537 


512 


509 


484 


E(7) 


594 


573 


567 


551 


Ax (4) 


598 


579 


575 


555 


E(8) 


700 


683 


694 


667 


B 2 (4) 


726 


701 


706 


667 


cell (see section IVD). This group of bands is slightly 
shifted downward with respect to experiments, as dis- 
cussed above. Our results compare fairly well with previ- 



ous lattice dynamics calculations of the DOS&^i and on 
dispersion relations.— However, sizable quantitative dif- 
ferences are present especially on angular dispersion with 
respect to shell model calculations^ To our knowledge, 
no experimental data are available on phonon dispersion 
relations. 



B. /3-Te0 2 

Phonons at the T-point can be classified according to 
the irreducible representations of the D 2 h point group as 
r = 9(A g + A u + B Xg + B 2g + B 3g) +8{B lu + B 2u + B 3u ). 
The three translation modes have been omitted. The g- 
modes are Raman active while B it-modes are IR active. 
Only B it-modes couple to the macroscopic longitudinal 
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field via the non-analytic part of the dynamical matrix. 
The calculated phonon frequencies at the T-point, ne- 
glecting the contribution of the longitudinal macroscopic 
field is given in Table HVBl The results refer to calcu- 
lations at the experimental lattice parameters with the 
BLYP functional which turned out to provide the best 
agreement with experiments for a-Te02 and for /3-Te02 
as well. 

The Raman tensor for the active modes of an or- 
thorhombic crystal has the formi^I 



A a 



B 2 (y) 



B 3 (x) 



f 



For /3-Te02 only the g-modes have to be considered. 

The calculated coefficients of the Raman tensor are 
given in Table [TVBl for all active modes. Experimental 
Raman spectra are available only for polycrystalline sam- 
ples in backscattering geometry for unpolarized light. 7 
The corresponding theoretical spectrum can be obtained 
by averaging the Raman cross section over the solid an- 
gle by making use of Eq. [7] The resulting theoretical 
Raman spectrum of /3-Te02 is compared with the exper- 
imental powder spectrum in Fig. 1131 The linewidth of 
the most intense Raman peaks in the theoretical spec- 
tra has been obtained by fitting the experimental spec- 
tra with a sum of Lorentzian functions. The results are 
given in Table IIV Bl The linewidth of the weaker or less 
identifiable peaks is set to 4 cm -1 . Overall the exper- 
imental spectrum is well reproduced, but for a redshift 
of the high frequency modes. The underestimation of 
the frequency of the stretching modes is probably due to 
the overestimation of the Te — O bond length (cfr. sec- 
tion IIIB). The frequency redshift with respect to experi- 
ments is even larger when the PBE functional is used, as 
occurs for the a-TeC>2 (cfr. Table HX|) . The theoretical 
Raman spectrum allows to assign the character of the 
modes responsible for the experimental peaks. More on 
the displacement pattern of the different modes is given 
in section IV D where the Raman spectra of the three 
phases are compared. 

The B it-modes are subject to angular dispersion due 
to the coupling with the macroscopic longitudinal field. 
The dependence of the phonon frequencies on the polar 
angles of the q vector approaching the T-point is given in 
Fig. HH IR oscillator strength of the active modes is re- 
ported in Tabic HVBl To our knowledge no experimental 
IR data are available for single crystals. To compare with 
experimental data available we have thus computed the 
IR absorption spectrum for a polycrystalline sample as 
described in the previous section. The result is reported 
in Fig. 1121 it has to be compared with the experimental 



FIG. 13: Unpolarized Raman spectra of polycrystalline a-, /3-, 
and 7-Te02 for backscattering geometry. Theoretical spectra 
are calculated with the BLYP functional at the experimen- 
tal lattice parameters. Experimental spectra are also shown 
for comparison.- For each phase the experimental spectrum 
is normalized over the strongest high-frequency peak of the 
theoretical spectrum. 
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spectrum reported in Fig. 2 of Rcf. 8 . The experimental 
peaks are reproduced qualitatively, but for an inversion 
in the relative intensities of the two peaks around 600 and 
700 cm -1 . Moreover, the theoretical spectrum shows a 
shoulder at high frequency which is absent in the exper- 
imental spectrum. Our choice of the phonon linewidth 
might be partially responsible for these disagreements 
with the experimental data. In fact, in the lack of ex- 
perimental data, the linewidth of the phononic modes 
is chosen similarly to those of a-Te02, i.e. 5 cm -1 for 
frequency below 150 cm -1 and 15 cm -1 for higher fre- 
quencies. 
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TABLE X: Theoretical phonon frequencies of (3-Te<D 2 at the T point, oscillator strengths (fj in Eq. [4} of IR active modes and 
squared coefficients of the Raman tensor of the Raman active modes (in units of 10~ 3 A 3 , see section IV). The contribution of 
the inner longitudinal macroscopic field is not included (LO-TO splitting) . All the theoretical values correspond to calculations 
with the BLYP functional at the experimental lattice parameters. The phonon linewidths (Aw) are obtained by fitting the 
experimental Raman spectra of a poly crystalline sample^ with a sum of Lorentzian functions (Fig. 1 1 3f) . 



Mode 


u) (cm 1 ) 


fi 


a 2 (d 2 ) 


b 2 (e 2 ) 


c 2 (f) 


Alu 


Mode 


lo (cm 1 ) 


fi 


a 2 (d 2 ) 


b 2 (e 2 ) 


c 2 (f) 


A u (l) 


38 












B 3u (5) 


310 


2.893 








B 3u {l) 


44 


0.049 










Bi g (5) 


312 




0.171 






Bi«(l) 


55 


0.092 










B lu (5) 


313 


6.503 








B 39 (l) 


56 








3.828 




B 39 (5) 


317 








0.002 




67 




0.576 


1.607 


2.645 




B 2u (5) 


325 


0.940 








Bi g (l) 


72 




0.000 








A u (6) 


336 










B 2g W 


77 






0.310 






#3 9 (6) 


349 








0.000 


A u {2) 


85 












Bi 9 (6) 


352 




0.400 






B 2u {l) 


89 


5.487 










A g (6) 


354 




2.525 


0.078 


0.120 


M2) 


108 




0.542 


3.133 


1.998 




Au(7) 


377 










B 2g (2) 


123 






0.186 






B 2u (6) 


390 


0.455 








B lu {2) 


144 


5.503 










B 2g (6) 


397 






0.127 




B 3u (2) 


153 


0.586 










A g (7) 


404 




0.112 


1.411 


4.031 


Bi a (2) 


159 




0.612 








B 2g (7) 


408 






0.129 




Sx„(3) 


161 


0.398 










Bin (6) 


421 


0.590 








A u (3) 


163 












B 3u (6) 


424 


0.000 








B 3g {2) 


169 








2.952 


3.5 


Bi 9 (7) 


452 




0.020 






B 2 u(2) 


170 


0.015 










B Sg (7) 


457 








0.306 


B 3u (3) 


171 


0.478 










B 3 u(7) 


518 


0.115 








4,(3) 


180 




1.157 


2.548 


1.650 


3.2 


B lu (7) 


520 


1.776 








B 2g (S) 


188 






0.324 






A g (8) 


535 




0.579 


1.192 


5.887 


B 39 (3) 


208 








0.044 




B 2g (8) 


541 






0.536 




Bi s (3) 


209 




0.003 








A u (8) 


546 










4(4) 


211 




1.163 


3.532 


2.730 


3.4 


B 2u (7) 


548 


0.000 








B 29 (4) 


212 






0.257 






B 3g (8) 


591 








8.965 


B 29 (5) 


222 






0.049 






Bi g (8) 


593 




1.586 






MS) 


223 




7.267 


0.819 


1.616 


4.4 


A g (9) 


614 




6.423 


18.540 


13.575 


Bi 9 (4) 


231 




0.138 








B 2g (9) 


615 






0.724 




Ba 9 (4) 


234 








0.001 




B 3u (8) 


640 


0.618 








A»(4) 


235 












A u (9) 


668 










B 2u (3) 


247 


0.666 










B 2u (8) 


670 


1.277 








B 2u (i) 


254 


3.949 










B 3g (9) 


675 








0.475 


A u (5) 


258 












Bi g (9) 


675 




0.023 






B 3u {<i) 


265 


0.094 










Bi«(8) 


706 


0.000 








B lu (4) 


267 


2.050 























C. 7-Te0 2 

Phonons at the T-point can be classified according to 
the irreducible representations of the D2 point group as 
r = 8(B 1 + B-2 + B 3) ) + 9 A. The three translational 
modes have been omitted. The calculated phonon fre- 
quencies at the T-point, neglecting the contribution of 
the longitudinal macroscopic field is given in Table 1151 
The results refer to calculations at the experimental lat- 



tice parameters with the BLYP functional. All the modes 
are Raman active while only B-modes are IR active and 
display angular dispersion due to the coupling with the 
longitudinal macroscopic field (Fig. [T5|) . IR oscillator 
strength of the active modes is reported in Table IXll Still 
no experimental IR data are available for 7-TeC>2 single 
crystals. To compare with experimental data available 
we have thus computed the IR absorption spectrum for 
a polycrystalline sample as described in the previous sec- 
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tion. The result is reported in in Fig. [T^J it has to be 
compared with the experimental spectrum reported in 
Fig. 2 of Ref A In the lack of experimental data, the 
linewidth of the phononic modes is chosen similarly to 
5 cm -1 for frequency below 150 cm -1 and 15 cm' 1 for 
higher frequencies, analogously to a-Te02- The theo- 
retical spectrum reproduces all the main features of the 
experimental spectrum^, but for a redshift of the theo- 
retical modes at high frequency as occurs for a-Te02. 
The Raman tensor for active modes has been given in 
the previous section. 

The coefficients of the theoretical Raman tensor are 
given in Table XI for all the active modes. As for (3- 
TeC>2, experimental Raman spectra are available only for 
polycrystalline sample in backscattering geometry. The 
corresponding theoretical spectrum can be obtained by 
averaging the Raman cross section over the solid angle. 
However, as opposed to the case of /3-TeC>2, Raman ac- 
tive modes in 7~Te02 can be subject to angular disper- 
sion (cfr. Fig. [15] ). As a consequence, formula Eq. 
[7] can not be used and the integral over the solid an- 
gle has to be performed over discrete angular values of 
the phononic wavevector q approaching the T point. For 
each wavevector q transfered in the scattering process 
the phonon eigenvector and frequency is computed by 
including the non-analytic part of the dynamical matrix. 



The solid angle average is performed over 256 angles in- 
dependent by symmetry. The calculated powder Raman 
spectrum of 7-Te02 is compared with the experimental 
spectra and with the spectra of the other phases in Fig. 
[T3l The linewidth of the theoretical peaks are fitted to ex- 
periments as discussed for /3-TeC>2 in the previous section. 
The linewidth of the weaker or less identifiable peaks is 
set to 4 cm -1 . The linewidth of the Raman active modes 
are reported in Table XI. The angular dispersion has a 
minor effect on the Raman spectrum since the most in- 
tense peaks turn out to have A character for which no 
coupling exists with the longitudinal macroscopic field: 
the Raman powder spectra of 7-Te02 including or ne- 
glecting the angular dispersion are indistinguishable on 
the scale of Fig. [13J The agreement between theory and 
experiments is overall good as for the other phases. Still 
a sizable redshift of the high frequency modes is present. 
Analogously to the other phases, the error is still slightly 
larger by making use of the PBE functional. Moreover, 
a very weak experimental peak around 790 cm -1 has no 
counterpart in the calculated spectrum. Although three 
Raman active modes fall in that region, none of them has 
an appreciable Raman intensity. More on the displace- 
ment pattern of the different modes is given in section 
IV D where the Raman spectra of the three phases are 
compared. 
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TABLE XI: Theoretical phonon frequencies of 7-TeC>2 at the 
r point, oscillator strengths (fj in Eq. U) of IR active modes 
and squared coefficients of the Raman tensor of the Raman 
active modes (in units of 10 -3 A 3 , see section IV). The con- 
tribution of the inner longitudinal macroscopic field is not 
included (LO-TO splitting). All the theoretical values cor- 
respond to calculations with the BLYP functional at the ex- 
perimental lattice parameters. The phonon linewidths (Aw) 
are obtained by fitting the experimental Raman spectra of a 
polycrystalline sample 7 - with a sum of Lorentzian functions 

(Fig. m- 



Mode 


lo (cm 1 ) 


fj 


a 2 (d 2 ) 


b 2 (e 2 ) 


c 2 if) 


Auj 


A(l) 


52 




0.029 


0.813 


0.137 




B 2 (l) 


90 


8.215 




0.002 






S 3 (l) 


111 


0.549 






0.084 




A(2) 


122 




2.222 


1.502 


0.025 


4.2 


Bi(l) 


135 


0.307 


1.190 






3.7 


Si (2) 


137 


1.627 


1.242 






3.7 


B 2 (2) 


142 


0.069 




0.076 






A(3) 


158 




0.693 


0.612 


2.512 


4.5 


B 3 (2) 


164 


0.144 






0.000 




A(4) 


169 




1.138 


2.259 


2.875 


4.5 


B 3 (3) 


183 


5.049 






0.003 




S 2 (3) 


200 


0.141 




0.117 






Si (3) 


201 


1.362 


0.017 








A(5) 


212 




0.311 


0.455 


0.112 




S 3 (4) 


265 


0.517 






0.264 




S 2 (4) 


275 


1.935 




0.323 






Si (4) 


280 


2.097 


0.457 








A(6) 


291 




0.480 


0.529 


0.004 




S 3 (5) 


293 


2.036 






0.661 




Bi(5) 


305 


0.037 


0.387 








S 2 (5) 


353 


0.608 




0.040 






A(7) 


410 




1.467 


0.794 


29.271 


13.4 


Si (6) 


414 


1.670 


0.601 








S 2 (6) 


415 


3.866 




0.257 






B 3 (6) 


496 


0.895 






0.119 




A(8) 


575 




0.033 


0.151 


9.968 


12.4 


Bi(7) 


585 


1.519 


0.231 








B 2 (7) 


633 


1.147 




0.006 






A(9) 


636 




10.494 


16.807 


3.977 


8.7 


B 3 (7) 


657 


0.577 






0.036 




B 3 (8) 


699 


0.158 






0.026 




B 2 (8) 


737 


0.088 




0.007 






Si (8) 


768 


0.005 


1.955 






4.5 



D. Comparison of Raman spectra of the three 
phases 

In this section, we discuss the assignment of the ex- 
perimental Raman peaks to specific phonons and com- 
pare the displacements patterns of these phonons in the 



three different phases. Firstly, we must say that the ab- 
initio results presented here are overall in good agree- 
ment with previous results from lattice dynamics calcu- 
lations with empirical interatomic potentials'^ ' 10 ' 11 The 
displacement patterns of the phonons responsible for the 
main Raman peaks of a-, /?-, and 7-Te02 are reported in 

Figs, m mi na 

As pointed out in Refiii the T-point phonons of a- 
Te02 can be interpreted in terms of the vibrations of 
the Te02 molecule given in Table I. The modes Ai(4) 
and £>i(5) responsible for the strong Raman peak at 598 
cm -1 and the weak peak at 537 cm -1 (cfr. Fig. [T3]) corre- 
spond to symmetric and antisymmetric stretching modes 
of the Te02 molecular units (Fig. [T6"b-d). The mode 
Ai(3) responsible for the Raman peak at 383 cm -1 (Fig. 
IT3"| is a bending mode of the TeC>2 molecule (Fig. [TBb). 
The mode E(6) responsible for the weak Raman peak at 
336 cm" 1 is a librational mode of TeC>2 molecule which 
modulates the length of the two axial bonds (Fig. n"6"a). 
The Raman peaks at lower frequencies (120-150 cm" 1 ) 
are due to librational modes of the TeC>4 units (-E'(l), 
Ai(l), i?2(l) in Table VIII). The displacement pattern 
shown in Fig. [T7J] are in good agreement with those of 
the corresponding modes obtained in Refi 7 ^ from lattice 
dynamics calculations with empirical interatomic poten- 
tials. However, the relative Raman intensities are bet- 
ter reproduced by our ab-initio calculations than by the 
Bond-Polarizability Model (BPM) used in Ref£±i. 

Concerning /3-Te02, the mode B 3g (8) responsible for 
the Raman peak at 591 cm" 1 is an antisymmetric stretch- 
ing mode of the TeC>2 molecular unit (Fig. [T7b, in 
agreement with Refi^i). The mode A g (9) responsible 
for the Raman peak at 614 cm" 1 is instead a stretch- 
ing mode of just one (the shorter) of the two equatorial 
bonds (Fig. [T7p!). The mode A g (8) at 535 cm" 1 is the 
stretching mode of the other longer equatorial bond (Fig. 
ITTb). This is in contrast with the results of lattice dy- 
namics calculations with empirical potentials of Refi^i, 
which gives a mixing of the two modes to produce a sym- 
metric and antisymmetric stretching mode of the TeC>2 
molecule. The mode A g (7) at 404 cm -1 is a vibration of 
the Te2C>2 ring which modulates the length of the axial 
bonds, in agreement with Refi^i. However, the bond- 
polarizability model used in Ref. 7 ' 11 gives unsatisfactory 
results for the relative Raman intensities. 

Concerning 7-Te02, the Raman peaks at 768 and 636 
cm" 1 correspond to stretching modes of the shortest 
Te-0 bond in-phase (A(9)) or out-of-phase (i?i(8)) 
within the four TeC>2 molecules in the unit cell, respec- 
tively (Fig. [THb-d). The modes A(8) and A(7) responsi- 
ble for the Raman peaks at 575 and 410 cm" 1 are mixed 
stretching and bending modes of the Tc — 02' — Te' 
bridge with the Te atoms fixed (Fig. [T8"a-b). Again, the 
displacement pattern of these modes are in good agree- 
ment with those calculated in Ref'^ from empirical 
interatomic potentials. The analysis of the vibrational 
spectra of a-Te02 and /3-Te02, presented above, identify 
clearly the TeC>2 molecules as the building units of the 
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FIG. 16: (Color online) Displacement pattern of the phonons responsible for the main Raman peaks of polycrystalline a-TeC>2 
(cfr. Fig. I13|l . The TeC>2 molecules formed by the two shorter equatorial bonds are shadowed. 




crystal. However, the interaction between the molecules 
is large. For instance, the symmetric stretching modes 
of the Te0 2 unit shifts from 598 cm" 1 (Al(4)) to 726 
cm -1 (£? 2 (4)) depending on the phase relation between 
the motion of the different molecules. As pointed out 
in Refiiii, the 7-TeC>2 phase behaves differently with re- 
spect to a- and /3-TeC>2 for what concerns the vibrational 
properties. Indeed, the displacement patterns reported 
in Fig. HU suggest a chain-like structure, the most in- 
tense Raman peaks consisting of stretching modes of the 
Te - 02' - Te' bridge and of the shortest Te - O bond 
which might be seen as side group of the polymeric chain. 
Still, the chains are strongly interacting since the A(9) 
and -Bi(8) modes which consists of the same intrachain 



vibration with different interchain phase relation differ in 
frequency as much as 130 cm" 1 . We note that in the j3 
phase the lowest energy phonons (A u (l), B 3u (l), Bi u (l) 
at 38-55 cm" 1 ) correspond to rigid out of-phase transla- 
tional modes of the two layers in the unit cell. In 7-Te02, 
the lowest mode at 52 cm" 1 corresponds to a rigid trans- 
lation of the chains along the chain axes, out-of- phase 
between the two chains in the unit cell. The other two 
modes at higher frequency (80, 111 cm -1 ) involve instead 
sizable deformations of the chains themselves. In a-TeC>2 
the low frequency modes (below 140 cm" 1 ) can be well 
described as rigid motions of the TeC>4 units around the 
bridging oxygen atoms. 
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FIG. 17: (Color online) Displacement pattern of the phonons responsible for the main Raman peaks of polycrystalline /3-Te02 
(cfr. Fig. [13}. The TeC>2 molecules formed by the two shorter equatorial bonds are shadowed. 




C) 8^(8), 591 crri 1 d) \{9), 614 cm 1 
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FIG. 18: (Color online) Displacement pattern of the phonons responsible for the main Raman peaks of polycrystalline 7-Te02 
(cfr. Fig. [13}. The Te02 molecule formed by the two shorter equatorial bonds are shadowed. 




a) 





A(7),410cm 1 b) 



A(8), 575 cm 



-1 




c) 



A(9) 1 636 cm 



d) 



B^S), 768 cm 1 
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V. CONCLUSIONS 

Based on first principles calculations, we have stud- 
ied the structural and vibrational properties of the three 
crystalline phases of TeC>2 {a, (3, 7). Phonon dispersion 
relations and IR and Raman spectra have been computed 
within density functional perturbation theory. The cal- 
culated Raman and IR spectra are in good agreement 
with available experimental data, and with previous lat- 
tice dynamics calculations based on empirical interatomic 
potentials ^ 10 ! 11 , but for an underestimation of the fre- 
quency of the Te-0 stretching modes (above 500 cm -1 ). 
The modes at high frequency (400-800 cm -1 ) can be well 
described in terms of vibrations of TeC>2 molecular units 
in a- and /3-Te02 which would then be seen as molecu- 
lar crystals. The analysis of the vibrational spectrum of 
7-TeC>2 suggest instead a stronger intermolecular inter- 
action supporting the picture of 7~Te02 as a chain-like 



structure, more connected than a- and /3-TeC-2. How- 
ever, the calculated elastic constants and electronic band 
structure do not show strong anisotropies which would 
have been expected from the description of 7-TeC>2 as the 
assembling of aligned polymers proposed experimentally. 
As a matter of fact, TeC>2 molecules in a- and /3-Te02 and 
polymerized TeC>2 chains in 7-Te02 are strongly interact- 
ing, possibly via electrostatic coupling, as suggested by 
the strong ionic character of the Te-0 bond emerged from 
previous ab-initio calculations . 21 ' 25 
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